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We propose a new mechanism for generating a single transverse spin asymmetry (STSA) in polar- 
ized proton-proton and proton-nucleus collisions in the high-energy scattering approximation. In 
this framework the STSA originates from the q — > qG splitting in the projectile (proton) light-cone 
wave function followed by a perturbative (C-odd) odderon interaction, together with a C-even inter- 
action, between the projectile and the target. We show that some aspects of the obtained expression 
for the STSA of the produced quarks are in qualitative agreement with experiment: STSA decreases 
with decreasing projectile xf and is a non-monotonic function of the transverse momentum kr- In 
1 our framework the STSA peaks at kr near the saturation scale Q s ■ Our mechanism predicts that the 

quark STSA in proton-nucleus collisions should be much smaller than in proton-proton collisions. 
f"^ ■ We also observe that in our formalism the STSA for prompt photons is zero. 
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I. INTRODUCTION 



In recent years a significant theoretical effort has been directed toward understanding the origin and the interesting 
properties of transverse spin asymmetries observed in high-energy scattering experiments. The single transverse spin 
asymmetry (STSA) in polarized scattering P'(jp) + A — > h(k) + X is an observable describing the correlation between 
1^ the transverse spin vector S of the projectile P and the transverse momentum k of the produced hadron h. This 
[ correlation can be expressed either in terms of the asymmetry in the scattering of the spin-up and spin-down transverse 
spin states, or in terms of the "left-right" (fc O — fe) asymmetric momentum distribution in spin-up scattering: 



a 2 kdy y > d 2 kdy^ ' u\£-±u j 



a / t \ a^-K. ay a^K, ay a*K ay *■ > a*K ay v ' »^ j . . 

• Interpreted as a left-right asymmetry in the produced hadron distribution, the sign of An is fixed to reflect a pref- 
\ erential scattering of particles to the beam- left when An is positive and a preferential scattering to the beam-right 
■ when An is negative. More concretely, if the polarized projectile moves along the +z axis and the transverse spin 
is oriented along the +x axis, then positive An corresponds to more outgoing particles produced along the —y axis 
than along the +y axis. The geometry of polarized scattering is illustrated in Fig. [TJ The STSA can also be thought 
of as a spin- momentum coupling term proportional to (S x p) ■ k in the particle production cross section, where p is 
the momentum of the incoming projectile. For a review of the current status of STSA physics see [l[ . 
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FIG. 1. A sketch of the geometry of polarized scattering: the incoming transversely-polarized projectile is moving along the 
+2-axis, while its spin is pointing along the +x'-axis. The positive j/-axis points to the right of the beam. 

In the 1990's, the E581 and E704 collaborations at Fermilab reported data on hadron production in transversely 
polarized proton-proton scattering that showed large unanticipated transverse spin asymmetries of up to 30-40% 0- 
|7|. Qualitatively, the asymmetries were consistent with zero for mid rapidities, but increased rapidly in the forward 
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A N vs ICj, for 71 Production 

(STAR 2008) 
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• x F >0.4 
■ x p <-0.4 



FIG. 2. Experimental data on the pion single transverse spin asymmetry An as a function of Feynman-a; reported by E581 and 
E704 collaborations (graphically reconstructed from Q|, shown in the left panel) for 0.7 < kr < 2.0 GeV/c, and as a function 
of the pion transverse momentum kr collected by the STAR collaboration [l^] (right panel) . 



scattering direction, as illustrated by the data shown in the left panel of Fig. [2] More recently, the PHENIX, STAR, 
and BRAHMS collaborations at the Relativistic Heavy Ion Collider (RHIC) have studied transverse spin asymmetries 
over a wide kinematic range 

at = 200 GeV The data the y have presented confirmed and extended 

the Fermilab results, and also indicated a non-monotonic dependence of STSA on the transverse momentum fcx of 
the produced hadron [l2|, HH ■> shown in the right panel of Fig. [2j 

At the time these spin asymmetries were first observed, there was no theoretical framework to understand them; 
on the contrary, prevailing wisdom expected spin-dependent effects to become negligible at high energies fl4| . Much 
theoretical progress in classifying and modeling sources of STSA has been made since then, predominantly within the 
framework of collinear factorization. 

There are three stages in the scattering process where an asymmetry could be generated, illustrated in Fig. EI Sivcrs 
, [l6( proposed the existence of a correlation between the spin of the polarized hadron and the transverse momentum 
distribution of its partons. This Sivcrs effect would act within the parton distribution function (PDF) (or, equivalcntly, 
the initial state wave function) of the polarized projectile to generate an asymmetry. Being a part of the PDF, this is 
an intrinsically non-perturbative process. The asymmetry could also be generated in the partonic scattering processes 
themselves; since the leading-order scattering process does not gen erate an asymmetry, any contribution would come 
from "higher-twist" terms in the interaction, as suggested in [17H261 | . Such an asymmetry generated by the interaction 
may be perturbative, depending on the kinematics. Finally, Collins [27j proposed a similar correlation between spin 
and transverse momentum that could occur during hadronization, as a produced parton undergoes fragmentation into 
final-state hadrons. The Collins effect would couple the spin of a produced parton to the momentum dependence of 
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FIG. 3. Sketch of the three potential sources of asymmetry: in the parton distribution function via the Sivers effect, via partonic 
interaction, or in the fragmentation function via the Collins effect. 
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its fragmentation function, resulting in a left-right asymmetry within the final-state jet. Again, the Collins effect is a 
low-energy, intrinsically non-perturbative process. 

Recent studies from STAR [28| were able to resolve the angular distribution within the jets and thus probe the 
Collins function directly. Their data suggests that the Collins effect's contribution to the asymmetry is small or 
consistent with zero, as illustrated in Fig. H) 

Collins Effect: A N vs cos y in Jl° Production 

(STAR Preliminary 2011) 



• x F >0.3 
■ x p <-0.3 



<cos y > 



FIG. 4. Recent preliminary data from the STAR collaboration (graphically reconstructed from [2Sj) for the 7r° asymmetry 
within final-state jets as a function of the angle 7 between the outgoing pion and the jet thrust axis. The Collins contribution 
is proportional to the slope of the data, and is consistent with zero. 



We are interested in understanding which pcrturbative mechanisms within the interaction can give rise to transverse 
spin asymmetries, and we would like to study the interaction with proton-nucleus (jr + A) collisions in mind. When 
the target is a heavy nucleus at high energies, or when the energy of the proton-proton collision is high enough, 
the projectile scatters off of the target's small- a; parton distribution where parton densities are high. The effective 
interactions in this small- a: regime are well described by the parton saturation/ Color Glass Condensate (CGC) formal- 
ism [Il],[i|,[3l|,[3l-^ which resums the multiple rescatterings in these dense 
color fields and incorporates the small-x evolution of those fields. This saturation framework is expressed naturally in 
coordinate space and introduces a characteristic saturation scale Q s describing the color-charge density fluctuations 
within the target. 

With the CGC interactions that mediate p^ + A scattering in mind, it is natural to search for the asymmetry using 
the language of the small- a; interactions, light-cone perturbation theory (LCPT) Ell, instead of the usual collinear 
factorization. Using LCPT, we calculate spin-dependent processes that contribute to the projectile's light-cone wave 
function, then convolute this wave function (squared) with the known scattering amplitudes from CGC that are 
responsible for the interaction of partons with the target. Surprisingly enough this procedure yields a non-zero STSA 
for quark and gluon production in p^ + A scattering. Our approach yields a different but complementary picture to 
the usual collinear factorization mechanisms for generating STSA and may provide new insight into the underlying 
physics for generating transverse spin asymmetry. 

In the past the interplay of spin and small- a; evolution was investigated in [59| using the standard Feynman diagram 
approach. In more recent years the role of spin has also been investigated within the McLerran-Venugopalan (MV) 
model j3oT - l37| in [6(| and within Mueller's dipole model [32|-[34j in [6l|]. The transverse spin asymmetry in the CGC 
formalism was studied in [62| - 65] , with the STSA being a completely non-perturbative effect on top of the perturbative 
CGC dynamics. Other non-perturbative approaches include using QCD instantons to generate STSA as discussed in 

M- 

In this paper, we present an analysis of STSA generated by the interaction of the spin-dependent LCPT wave 
function of the projectile with the target gluon field in the saturation/CGC framework. The paper is organized 
as follows. In Sec. |TT] we derive the general expressions for the quark, gluon, and photon single transverse spin 
asymmetries in our formalism. These general results are given in Eqs. (|32|) . (|40|) and (|44|) . We show that to generate 
the asymmetry the projec tile needs to interact with the target via a C-odd scattering amplitude, commonly known as 
the QCD odderon |67H74| . In Sec. lIIIIw e simplify the interaction of the light-cone wave function with the target using 
the Glauber-Mueller approximation [75J , and we present a closed-form estimate of the asymmetry for produced quarks. 
We find that a C-odd exchange alone is insufficient to generate the asymmetry, and needs to be accompanied by a 
C-even exchange for the asymmetry to be non-zero. We also see that the generated STSA (An) is a non-monotonic 
function of transverse momentum kx and an increasing function of xf, in qualitative agreement with the data shown 
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in Fig.[2j Our rough estimates of the asymmetry are plotted in Figs. [T3l[T4l and [T5l which, when compared to the data 
in Fig. [3J suggest that it may even be possible to achieve quantitative agreement with the data using our approach in 
a detailed phenomcnological study. We conclude by summarizing our main results in Sec. IIVI 



II. GENERAL RESULT: COUPLING SPIN TO INTERACTION C-PARITY 



In this paper we consider two processes contributing to the STSA of hadron production in the scattering of trans- 
versely polarized protons on an unpolarized target (proton or nucleus), p^ + A — > h + X: quark production and gluon 
production. We will also consider STSA in prompt photon production, p^ + A — > 7 + X. The proton's transverse 
spin is transmitted to a quark in its wave function (giving rise to the transversity distribution); for simplicity we 
consider the scattering of this transversely polarized quark on the target separately from other (spectator) quarks in 
the proton. This process is illustrated in Fig. [5j where the high energy interaction between the projectile quark and 
the target nucleus is schematically denoted by gluon exchanges. 




nucleus, A 



q, G, 7 



> x 



J 



FIG. 5. Transversely polarized quark scattering in the field of the nucleus producing either a quark (q), gluon (G), or a prompt 
photon (7) along with extra hadrons denoted by A: + A — > (q, G, 7) + X. 



To introduce the methodology, we will concentrate on the quark production process, + A — > q + X. We 
are interested in isolating the spin-dependent contribution of the quark-target scattering. Certainly keeping only 
the eikonal interaction would be insufficient to reach this goal, since the eikonal scattering is independent of the 
quark polarization. A non-eikonal correction has to be included: in the multiple-rescattering Glauber-Mueller (75j 
approximation, the non-eikonal rescattering corrections are suppressed by powers of energy and are very small. A 
much larger spin-dependent contribution comes from the non-eikonal splitting of the projectile quark into a quark and 
a gluon, q — > qG, which is suppressed only by a power of the strong coupling a s . In the LCPT language the q — > qG 
splitting may take place either before or after the interaction with the target, as shown in Fig. [51 (Splitting during the 
interaction with the target is suppressed by powers of energy [76j].) The lowest-order diagrams shown in Fig. [5] that 
contribute to STSA in q^ + A — > q + X contain the emission of a single gluon from the polarized quark, where both the 
gluon and quark can scatter in the field of the target. (Multi-gluon non-eikonal emissions are also possible, but they 
are higher-order in a s and, hence, outside of the leading-order precision of this work.) The spin dependence of the 
process illustrated in Fig. [S] originates within the light-cone wave function of the quark-gluon system, which couples 
to the interaction in a way that generates the asymmetry. In this Section, we will first outline the calculation of the 
q — > qG light-cone wave function using LCPT. Then we will combine the resulting splitting wave function squared 
with the quark and gluon interactions in the field of the target and identify the contribution to the asymmetry. In the 
end we obtain general expressions for quark, gluon, and photon STSA's in our formalism. The incoming light quark 
has a particular flavor /; multiple quark flavors can be incorporated into our formalism by convoluting the obtained 
cross sections with quark distributions corresponding to different flavors (inserting the appropriate quark masses into 
our results below). 

Throughout this paper, we will work in light-cone coordinates p^ = (p + ,p~ , p) with p = p°±p 3 and the transverse- 
plane vector p = (p 1 ,p 2 ). Note that p ■ q = (1/2) (p + q~ + p~q + ) p ■ q and p^ p^ L = p + p~ — p\ with pt = We 
assume that the incoming projectile quark is moving along the light-cone x + -axis, while the target is moving along 
the x~-axis. We will work in the light-cone gauge of the projectile, A + = 0. 
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FIG. 6. Two contributions to the amplitude for the high energy quark-target scattering in LCPT. 

A. Light-Cone Wave Function and Transverse Polarization 

Consider the splitting shown in Fig. [7] of a transversely polarized quark with momentum p and polarization x — ±1 
decaying into a gluon (with momentum p — k, polarization A, and color a) and a recoiling quark (with momentum 
k and polarization %'). The projectile quark is traveling along the light-cone ic + -direction and the recoiling quark 
carries a fraction 



k+ 
p+ 



(2) 



of the incoming quark's longitudinal momentum. We do not restrict ourselves to the case of an eikonal quark emitting 
a soft gluon (1 — a -C 1), but work in the general case when both the quark and the gluon can carry comparable 
longitudinal momenta. 



p — k, A 



P, X 
— > — 




FIG. 7. The light-cone wave function for the q — > qG splitting. Vertical dotted line denotes the intermediate state. 



The calculation of the light-cone wave function corresponding to the diagram in Fig. [7] is different from other similar 
calculations in the literature (see e.g. [z3, [zl|) only in that now the incoming quark is polarized transversely. To 
account for the transverse spin of the projectile traveling in the x 3 = z direction, we define the axis of spin quantization 
to be the x 1 = x axis, and we need to construct the Dirac spinors corresponding to the spin projection eigenvalues 
X = ±1 on the x-axis. The proper Lorentz-covariant operator describing the spin orientation is the Pauli-Lubanski 
vector 



W u 



(3) 



where 



gup 



are the generators of the Lorentz group and we are using the convention in which £0123 = +1- 

For our purposes, we need to find the eigenspinors of W\ for a particle moving along the z-axis. Diagonalizing W\ 
in terms of the pure helicity eigenspinors gives (cf. e.g. [79| ) 



1 



V2 



(+*) 



xu { 



(-*)} 



(4) 
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where U± z are the spinors in the hclicity basis and x 
of W\ and of the Dirac operator: 



±1. These particular spinors U x are simultaneous eigenstates 



WiU x 
(7 • p - m) U x 



m 







Since the incoming quark is polarized transversely, it is convenient to use the same spinor basis (j4|) for the outgoing 
quark in Fig. [7J as well. Using the standard rules of LCPT [H| an d working in the transverse polarization basis 
^ for the spinors, we evaluate the light-cone wave function shown in Fig. [7] as 



^\ xx 'ik,p,a) = — 



(p-k)- 



VAQ c * u x(p) 



where 



= 0, 



2e A • (p-fc ) 

+ -k+ 



(5) 



(6) 



is the gluon polarization vector with the transverse components e\ = (—1/^/2) (A, i). 

In arriving at Eq. ((SJ) we have used the fact that the incoming state in Fig.[7Jcontains only the quark with momentum 
p while the intermediate state contains the quark and the gluon, as denoted by the vertical dotted line in Fig. [7J For 
diagrams where the polarized quark scatters in the nucleus before the gluon emission, as shown in the right panel of 
Fig. [6l the roles are reversed: the quark line p is the intermediate state, and the quark-gluon system is the final state. 
Since ^ init P~[ = YlfinalPT) ^ ne ener gy denominator reverses sign for final-state splittings. 

Using the on-shcll conditions explicitly gives the terms entering the energy denominator: 



Pt 



(P - kY 



(P-fc) 2 

p+ - k+ 



(7) 



We use the spinors in the Brodsky-Lepage convention |57l.l5a|. which, for a particle moving alon g th e z-axis, become 
helicity eigenstates. The spinor matrix elements for the Brodsky-Lepage spinors are well known j57ll58|. Making the 
change of basis (|4|) gives the relevant spinor products as 



U x ,(k) , U x (p) 



u 



Ak)i U x {p) 



2 6, 



ap 
S 



+ 



(fc 



a Pa 



(1 



ap 1 



a) i m x 5 z2 
i e 1 -' (k J j_ — a p 3 ^) + (1 — a) m x S l1 



(8b) 



and the 7 matrix element does not contribute to 7 ■ e* x since e x = in the light-cone gauge. Here e 



12 



^21 



1, 

0. With the matrix elements ([8]) it is straightforward to evaluate the light-cone wave function ([5]) in 



momentum space, obtaining 

^\ xx >{k,p,<y) 



where 



(fe — ap) 2 + rh 2 
e* x -{k- ap) [{I + a) 5 XX > + A (1 - a) 5 X - X > 

m = (1 — a)m 



V2 



(1 - a)x(5 xx > -X6 Xt - x >) 



(9) 



is a natural effective mass parameter in the wave function and T a are the SU(iV c ) generators in the fundamental 
representation. 



1 Note again that U x (p) becomes a spinor for a transversely polarized particle only for p = 0: Eqs. (£Sj) give us the matrix elements for 
spinors related to the Brodsky-Lepage spinors via Eq. Q, which do not necessarily correspond to transverse polarizations in the general 
case. 
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Now we can Fourier transform the wave function to coordinate space 

d 2 k d 2 p . lk . (z _ x) ip .( x -u) I a 



ip Xxx ,(x,z,a;u) 



(2tt) 2 (2tt) 2 



(10) 



with the transverse coordinates defined in Fig. [7] Since the momentum-space wave function depends only on k — - ap, 
one of the two integrals can be performed to yield a delta function S 2 [(x — u) + a [z — as)]. Performing the remaining 
momentum integral in Eq. (|10[) yields 



i/) Xxx ,(x,z,a;u) 



g T a 

— — S 2 [(x — u) + a(z — x)] ir 

X (1 - ol) . _ 

Ko{m\z — x\ 



V2 



A<5, 



z - 


- X 


v- 


-x\ 


-x' 


}■ 



Kiifh \ z — x\) 



(l + a)5 XtX > +\(l-a)5 x - x > 



(11) 



It is useful to separate out the color factor T a and the delta function from the rest of the wave function (denoted by 
^f\ xx '), such that 



i> a x , (as, z, a; u) = T a 5 Z [(it - x ) - a {z - x)] Vx xx > (z - x, a) 



(12) 



Finally, we need to square the wave function and sum over the final particles' polarizations. Here we are interested in 
producing a quark with a fixed transverse momentum, while integrating over all transverse momenta of the produced 
gluon in Fig. |S1 According to the standard prescription [56l [7a |. for coordinate-space scattering amplitudes this 
means that the gluon's transverse coordinate x will be the same both in the amplitude and in the complex conjugate 
amplitude (since its transverse momentum is integrated over in the cross section), while the quarks have different 
transverse coordinates between the amplitude and the conjugate amplitude (since they are the observed particles). 
(See Fig. [TOl below for the illustration of the amplitude squared.) The "square" of the light cone wave function (fTTj) 
with the above rule for the quark and gluon transverse coordinates is illustrated in Fig. [5J 




FIG. 8. Light-cone wave function from Fig.[7]squared. The vertical straight line separates the wave function from its conjugate, 
while the cross denotes the quark that we tag on. The untagged gluon's coordinate x is unchanged, but the quark coordinates 
differ (z and u vs. y and w, as explained in the text.) 



The wave function (jllj) squared $ x as shown in Fig. [8] contains one contribution which is polarization-independent 
and another which is proportional to the quark polarization eigenvalue \ 

®x( z - x >y- x ' a ) = ^^xx'iz-X'®) *A XX '(y-^ a ) = ^unp{z-x,y-x,a) + x^ P oi{z-x,y-x,a) . (13) 

A,x'=±l 

Substituting the wave function (flT|) into Eq. (fL3|) and performing the sums gives the unpolarized part as 



2a 



s ~ 2 

m 



(1 + a 2 ) ( f T M"* I* - *l) Ki(m \y - x\) 

\z -x\\y- x\ 



+ (1 - a) 2 K (m \z - x\) K (m \y - x\) 



and the transversely-polarized part as 
2 a, 



$ P oi = m a (1 - a) 

7T 



K (fh \y - x\) Ki(m \z - x\ 



— 

\y - x\ 



K x {m \y - x\) K (m \z - x\ 



(14) 



(15) 
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Note that Q unp is a scalar under rotations in the transverse plane, whereas <& po i has an explicitly preferred x 2 -direction 
(i.e., the y-axis) since it "knows" about the transverse polarization of the incoming quark. The x 2 -axis can be written 
as the direction of the p x S vector, since the incoming quark with momentum p is moving along the z-axis, while 
being polarized along the x = x^axis, such that S\\ x\. We show that the unpolarized part of the wave function 
squared $„ p contributes to the unpolarized quark production cross section da unp , while the polarization-dependent 
part of the wave function squared $ po ; generates the spin-asymmetric cross section d(Aa). 



B. Spin, Asymmetry, and C-Parity in Quark Production 

Having computed the q — > qG light-cone wave function, we can now construct the scattering cross section by allowing 
the wave function to interact with the small- a; field of the target nucleus. It is well known [SI . l49l | that eikonal quark 
and gluon propagators in the background color field a can be correspondingly written as fundamental and adjoint 
path-ordered Wilson lines 



V X = V exp 



2 



dx + T a A~ a (x A 



= 0,x) 



(16a) 



U b x a =Vexp 



2 



dx + t c A- c (x + ,x~ = 0,x) 



(16b) 



where t a, s are the SU(iV c ) generators in the adjoint representation and the projectile is moving along the light-cone 
x + -axis. In essence, this means that the the projectile's transverse position is not altered during the scattering, and 
the effect of the target field is to perform a net SU(-/V C ) color rotation on the projectile. The Wilson lines resum these 
interactions and give the total phase of that color rotation. They are illustrated in Fig. [9l Note that the adjoint 
Wilson line U ba is real- valued. 
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FIG. 9. Wilson lines resumming scattering in the small- a; field of the target. The quark propagator is in the fundamental 
representation (top), and the gluon propagator is in the adjoint representation (bottom). 



The Wilson-line approach is quite generic: if the target gluon field is quasi-classical, as in the case of the McLerran- 
Venugopalan (MV) model |35l - [37| . then correlators of the Wilson lines resum powers of a 2 s A 1 / 3 corresponding to 
the Glauber-Mueller multiple-rescattering approximation (75j . Non-linear small- x evolution resumming powers of 
a s Y ~ a s In s can be included into the correlators of the Wilson lines through the Balitsky-Kovchegov (BK) (49i — 
[o3j and Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) jio| - |48j evolution equations. Thus 
expressing the interaction with the target in terms of the Wilson lines (fT6]) allows for several different levels of 
approximation for this interaction. 

The scattering amplitude for quark production is composed of two sub-processes: the splitting of Eq. (|11[) calculated 
in Sec. Ill Al and the Wilson line scattering (|16[) of the quark and the gluon in the field of the target. These elements give 
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FIG. 10. The cross section for quark production in the polarized quark-nucleus scattering. 



two distinct diagrams contributing to the scattering amplitude shown in Fig. [5] above. To find the quark production 
cross section we need to square the diagrams in Fig. |6l keeping the transverse momentum of the quark fixed, as 
depicted in Fig. [TO] As discussed above, this implies that the transverse coordinates of the quark are different on both 
sides of the cut. Just like in other similar calculations [76|, [zl] the q — > q G splitting may occur either before or after 
the interaction with the target, both in the amplitude and in the complex conjugate amplitude, resulting in the four 
different terms shown in Fig. 1101 

Using Fig. [TO] we can write down the expression for the color-averaged amplitude squared (M 2 ^ in terms of Wilson 
lines and the wave function responsible for the splitting, remembering to reverse the sign in the wave function for 
splitting occurring after the interaction: 



<M 2 > 



N, 



Tr 



E 

A,x' 



Tr 



UZ a Ul c + Tr 



(17) 



u. 



ho 



Tr 



n xx -v u ^ xx ,vi 



u. 



ab 



where N c is the number of colors, the traces are taken over the fundamental representation indices, and summation 
is implied over repeated adjoint color indices. Substituting Eq. (TTOj) into Eq. (fTT)) and using the identities 



jjba T a = y x T b y 



Tr \A T a B T a ] 



- TrA Ti B 
2 



1 

2N C 



■ Tr \A B] 



for arbitrary N c x N c matrices A, B, we find 

(A/ 2 ) = C F 5 2 [u - x - a (z - x)] S 2 [w - x - a (y - x)] <P x (z -x,y-x) l (q) 
where the factor responsible for the quark's interaction with the target, denoted by I^, is given by 

1 „ r„ „ + n „ r„ „ + n 1 



(18) 



(19) 



Tr [V s V2] Tr [V x v£\ 



1 



2N C C F 



Tr [V u V£] Tr [V*Vj] + 



2 iV 2 C F 
1 



Tr [V z Vi] 



2N}C F 



Tr [V u Vr\ 



(20) 



Here Cp = (iV 2 — 1)/2N C is the fundamental Casimir operator of SU(iV c ), and the angle brackets on the right denote 
averaging over the field configurations of the target. 
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Defining the S-matrix operator for a fundamental-representation color dipole by 

D wv = ±-Tr[V w V2] (21) 

we can rewrite 1^ more compactly as 

1^ = D z y + D u w — ^ D z x D x w + 2 q f ^ z w ~ 2 Cf ^ U X ^ X V 2 N Cf ^ / 

As we have already mentioned, this interaction with the target can be evaluated either in the Glauber-Mueller 
multiple-rescattcring approximation or using the JIMWLK evolution equation. 

The expression (f2"2"j) simplifies in 't Hooft's large- N c limit, in which the correlators of several single-trace operators 
factorize, such that, for instance, (D ux D xy ) = (D ux ) (D xy ) [HI HH, HU . Defining 

D xy = (b xy ) = ±- (Tr [V x Vj] ) (23) 



we rewrite Eq. (f2"2"|) in the large- N c limit as 

large-7V c 



D z y ~\~ Du w D z x D x w D u x D x y . (^4) 



To compute the quark production cross sections, we need to Fourier transform the coordinate space amplitude 
squared of Eq. ([!§]) back to momentum space and include the appropriate kinematic factors. This is accomplished by 



da^ 1 



<Pkdy q 2(2tt) 3 I- a 



d 2 x d 2 y d 2 z d 2 u d 2 w e - ik <*- y) e lp - {u - w) (M 2 ) (25) 



with k and y q the transverse momentum and rapidity of the produced quark. Integrating over the delta functions in 
(|19|) imposes the kinematic constraints 

u = x + a (z — x) (26a) 

w = x + a (y — x) (26b) 

which relate the quark coordinates before and after the q — > qG splitting and describe the non-eikonal quark recoil. 
To make the incoming quark transversely polarized we need to put the transverse momentum of the incoming quark 
to zero: p = 0. We thus obtain the general result for the quark production in the q^ + A scattering 



da^ C F 



d?kdy q 2(2tt) 3 1 



J d 2 xd 2 yd 2 ze- lk(z - y) <P x (z-x, y-x,a)l M {x, y , z) (27) 



with $ x from Eq. CE3) and from Eq. ||22J). The expression (|2"T)) contains multiple rescatterings and non-linear 
small-x evolution between the projectile and the target. Note that it does not resum the small- a; evolution between 
the produced quark and the projectile (which can be included following f78j|). and hence is not valid for very small 
a (i.e., the values of a are restricted by a s lnl/a <C 1). Since, as we will see below, both the experimental STSA 
and STSA resulting from our production mechanism fall off with decreasing a, the region of interest in this work 
corresponds to a not being very small, where Eq. (|2T[) is fully applicable. 

Now we are in a position to analyze the symmetry properties of the wave function and the interaction. There are 
two relevant symmetries to consider: k — > k or "/c^-parity" and \ ~^ ~X spin-flip. Under fcr-parity, the quark and 
antiquark coordinates z and y get interchanged, z y. (Note that, due to Eqs. (|26l) . this also implies that u <H- w.) 
From our previous calculation of the wave function (|14[) . (|15[) , we note that $„, lp is spin- independent, a scalar under 
rotations in the transverse plane, and fcr-even. Similarly, \ ^poi is odd under spin-flip, a vector under transverse 
rotations, and fcr-even. We can explicitly (anti-)symmetrize the interaction with the target under z -H- y interchange 

Zsymm / anU = 7J ( I ± (z ^ y) J • (28) 
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To do so wc first need to decompose each dipole S- matrix into the even and odd pieces under the exchange of transverse 
coordinates, which corresponds to the C-parity operation exchanging the quark and the anti-quark [Til. [72j: 

D xy = S xy + iO xy (29a) 

S x y=^{D xy +D yx ) (29b) 

X y — (D X y Dy X ^j . (29c) 

The C-even real part of the target field-averaged S'-matrix S xy = (S xy ) is responsible for the total unpolarized cross 
section of the dipole-target interactions. Its small- a; evolution is given by the BK/JIMWLK equations. The C-odd 
imaginary part of the target-averaged S'-matrix O xy = ( O x v ) is known as the odderon interaction [6714691 ]. The 
small- a; evolution equation for O xy was constructed in [7lU73l | . and, in the linear approximation, was found to be 
identical to the dipole BFKL equation [32[ with C-odd initial conditions. The intercept of the linearized odderon 
evolution was found to be ao — 1 = 0, in agreement with the solution of the Bartels-Kwiecinski-Praszalowicz (BKP) 
[1(3, equation for the odderon found in [70 [.For the current status of the experimental searches for the QCD 
odderon and for an overview of the theory see 69( . 

With these explicitly symmetrized elements, it is straightforward to construct the symmetric and antisymmetric 
parts of the interaction with the target (|2"2"|) for quark production: 

-^syinm = |^S^ $ u w — - — (j3 z x S x w — O z x O x ui^j + - — — S z w 

N r /- A A A \ 1 



(^S u x S x y O u x O x y^j o J\[ (J ^ U ^ / ' (30a 



AWa A A A I _ 



-^anti — * ( y + O u w 2 [O z x S xw + S zx O x ™) + 2N (j^ O z -u 

Nr /A A A A \ 1 



2 C j? v » / 2 N C f 

In the large- N c limit these expressions simplify to 



U X S X y + S U X X y J + U y ) . (30b) 



T(9) 

symrri 



anti 



— Szy H~ S uw S Z x S X w S ux S X y H~ Ozx Oxw H~ Oux O xy , (31a) 

large- jv c 

2 [C z y + ^ O z x S x w O u x S x y S z x O x w S u x O x y\ . (31b) 

large- w c 

Because the production cross section preserves azimuthal rotational symmetry, any physical observable must behave 
the same way under both fcr-parity and spin flip; that is, the only terms which are non-zero are those which are even 
under the combined operation hr <& (spin flip). Note that $> U np is even under the spin-flip, while x^poi is odd under 
spin-flip, and both are even under fc^-parity. Hence in the $ x X^ 9 ' product in Eq. (|27[) the terms which give non-zero 

contributions to the cross section arc & U npZsymm and x^pol^alu- Tne other terms $ unp l^ lti and x^poi^iylnm 
contain a contradiction between their vector structure and their fc^-parity, so they vanish identically: for instance, 
the contribution to the cross section coming from ^unpl^anti ^ s oa -d un dcr /c^-parity and is a function of k\ only, 
which is possible only if it is zero. Similarly, the contribution of x^poi^-sylnm is even under fcy-parity and is a vector 
in transverse space, which implies that it is also zero. Of the nonzero terms, ^unp^sylnm does not change its sign 

under spin-flip; it generates the symmetric part of the distribution da unp in Eq. ([T]). At the same time x^poil^anti 
does change sign under spin-flip; it generates the transverse spin asymmetry we are looking for. 

Knowing these symmetry properties, we can explicitly construct the spin-dependent and spin-averaged cross sections 
d(Aa) and da unp for quark production from their definitions (JTJ) obtaining 

d(Aa^) = ^-^-J d 2 xd 2 yd 2 ze- lk ^ % ol {z - x , y - x,a) 2$«(s,l/, z) (32a) 

d<r { u % = ^ J d 2 xd 2 yd 2 ze-* k < z -^ $ anp (z -x,y-x,a) l^ mm (x , y , z) , (32b) 

where the wave functions squared are given by Eqs. (fT4|). (fT5|) . and the interactions are given by Eqs. ([30]) (and by 
Eqs. ([31]) in the large- N c limit). Eqs. (|32|) are one of the main results of this work: together with Eq. ([lj they give the 
single-transverse spin asymmetry An generated in quark production by the C-odd CGC interactions with the target. 
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The mechanism for the generation of the STSA in Eqs. (|32|) is different from both the well-known Sivers [ID, 
and Collins Jl?} effects. It appears difficult (if not impossible) to absorb the interactions of Fig. [TU] into the projectile 
wave function (distribution function): hence our result is different from the Sivers effect. In the above calculation 
the asymmetry is generated before fragmentation; hence the STSA resulting from Eqs. (f3"2"j) cannot be due to Collins 
effect either. As we will see below, the non-zero part of Eq. (|32a[) stems from the multiple interactions with the target 
(higher-twist effects), and its contribution is in fact zero in the linearized (leading-twist) approximation. In this sense 
the above mechanism for generating STSA is similar in spirit to the higher- twist mechanisms of [l7l - l26j . though a 
detailed comparison of the diagrams appears to indicate that the two approaches are, in fact, different. 

We have shown explicitly that the single-transverse spin asymmetry An occurs in the CGC framework as a coupling 
between the transverse spin of the projectile and a C-odd interaction with the target, driven by the odderon0 Note 
that to date there is no unambiguous experimental evidence for the QCD odderon. If our mechanism for generating 
STSA can be isolated experimentally from other contributions, it may constitute the first direct observation of the 
QCD odderon! To make such a distinction possible, one needs to determine phenomenological characteristics of our 
mechanism, such as its rapidity, energy, and centrality dependence; some of this work will be carried out below, while 
the rest, along with a proper phenomenological implementation of our results, is left for future work. 

Within LCPT, the real- virtual cancellations or "crossing symmetry" embodied in Eq. (|18[) allow particles appearing 
in the complex conjugate amplitude M* to be rewritten as their charge-conjugate particles appearing in the amplitude 
M . This is what gives rise to the natural dipole degrees of freedom within CGC, and in terms of the asymmetry, this 
feature translates the fcy-parity of the cross section into the C-parity of the interaction with the target. This is the 
reason the odderon appears naturally in the expression for the asymmetry. 

In the literature it is often emphasized that the STSA An is odd under the time reversal transformation T [27j . 
To elucidate how the obtained result transforms under T we note that T reverses the directions of 5, p, and k: 

S—¥—S, p—^—p, fc4 — k. (33) 

Since (px S) \\ xi is invariant under T, the only effect of time-reversal on the spin-dependent cross section (|32a[) is 

d(A(T^)(fc) 4 d(Aa^)(-k), (34) 

since it is a real Lorentz-scalar momentum-space quantity. Therefore time reversal is equivalent to the fcr-parity 
transformation discussed above, and, hence, to the z -H> y interchange. Because of the odderon exchange, the 
interaction with the target X^ ti (x , y , z) is z <-> y anti-symmetric, and is, therefore, T-odd. This results in the 
spin-dependent cross section being T-odd too, 

d{Aa {q) )(k) 4 d{Aa {g) )(-k) = -d(A C r ( " ) )(fc), (35) 

leading to the T-odd STSA An, in agreement with the standard expectations (27j . 

It is interesting to note that in the high-energy approximation considered here the application of time reversal to 
dipole correlators is equivalent to the application of C-parity, such that the STSA arises from the odderon exchange, 
which is both T- and C-odd. As one can check explicitly, under time reversal Wilson lines transform as 

14 4 TV <B T- 1 = (V^)* = Vj (36) 

with T in the superscript denoting transposition. Note that if the original Wilson line was a future-pointing integral 
along the x + -axis, the time-reversed Wilson line on the right of Eq. (|31)f can be thought of as a complex conjugate of 
a past-pointing Wilson line or a transpose of the future-pointing Wilson line. The integration in either interpretation 
runs along the x~ axis: we will not show or refer to such x + <-> x" interchange explicitly, as it can be eliminated by 
a simple relabeling of the x 3 -axis. Since the averaging over the nuclear target is T-even, which can be inferred from 
[IS Hfil, the (coordinate-space) time reversal of the odderon amplitude O xy = (O xy ) applies only to the odderon 
operator, such that 

(O xy ) 4 -(6 xy ), (37) 

because (D xy ) is invariant under transposition and T is antilinear (i 4 — i). Thus we see that the T-odd odderon 
exchange leads to the T-odd STSA observable An- 

Finally, the reader may wonder whether the cross section in Eq. (|32ap is in fact non-zero. While it is very difficult 
to carry out the integration in Eq. Q32ap exactly, we instead will evaluate Eq. (|32ap approximately in Sec. IIII1 showing 
that the cross section and the corresponding STSA An are in fact non-zero. However, first we would like to derive 
the analogues of Eqs. (f3"2"| for gluon and prompt photon production. 



2 In the past, the relation between the odderon and the single and double transverse spin asymmetries was investigated in I82I486H in the 
pomeron and reggcon formalism. 
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C. STSA in Direct Quark Production 

Another channel for quark production is shown in Fig. [TT] and contains a "virtual" gluon correction, without the 
gluon being present in the final state. The gluon and quark in the loop may still interact with the target, as depicted 
in Fig. 1111 It can be shown that similar diagrams with the gluon loop located either completely before or after 
the interaction with the target do not contribute to the STSA and are, therefore, not considered here. The process 
illustrated in Fig. [IT] only leads to quark production at a = 1, and may seem to be negligible if, in order to avoid 
the projectile's fragmentation region, we constrain ourselves to a < 1 kinematics for the produced quark. However, 
remembering that for phenomenological applications one would have to convolute our quark production cross section 
Q32[) with the quark distribution in the projectile proton, we see that a = 1 production may still lead to non-negligible 
quark production at Feynman-cc less than one, possibly avoiding mixing with the proton's fragmentation region. 




c. c. 



FIG. 11. Direct quark production in the polarized quark-nucleus scattering. 

The calculation of the cross section for the process in Fig. [TTJ can be straightforwardly carried out along the lines 
specified above, yielding for p = 



da (q) direct (J 



d 2 kdy q (2 irf 



8(1 — a) / d xd yd zd ud vd we 



-ik-(z—y) 



1 

f da' , 

/ / r *v(i> — x, v — x, a ) 

J 47r(l-a') xy ' 



6 2 [(u -x)-a'(v- x)] 5 2 [{z - x) - a' (v - x)} 8 2 (y - w) l[ q) direc \ x j y j „) 



■ 5 2 [(w -x)-a'(v- x)} 8 2 [{y - x) - a' {v - x)} S 2 (z - u) dmc \ x jZ , v ) 



(38) 



with a 1 the longitudinal momentum fraction of the incoming quark carried by the quark in the loop, y q ~ In 1/a (up 
to an additive term defining the zero rapidity direction), and the interactions with the target 



X(q) direct / 
i (x,y,v) 



Nr b D 1 p ) 



r {q) direct 



(x, z, v) 



2C F 
N c 
2C F 



2N C C F 
1 

2N C C F 



(39a) 
(39b) 



To obtain the contribution to the numerator of STSA in Eq. (UJ from Eq. (|38p we keep the spin-dependent part 
of the wave function squared, $ po h an d anti-symmetrize the integrand with respect to the z o y interchange. This 
gives 



d(Aa {q) d%rect ) 



C F 



2(2tt) 3 



6(1 -a) / d 2 xd 2 yd 2 ze- lk{z - v) 



da' 



(1 - a') a' 



z — X z — X 



a' a' 



I \ t(q) direct / 
T > a X anU (x,y,Z)~ <$> p ol 



y X y X f \ (q) direct 

- a T„' 



a' a' 



(x,z,y) 



(40) 
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with the anti-symmctrized interaction given by 



f{q) direct 
anti 



{x,y,z)=I\ 
N, 



(q) direct 



x,y, x 



r (q) direct 



c. 



V ®x+-^ (z-a 



'-'a:, x+ (z— x) Ox, y 



N C C F 



x, y , x 

6 



y, x+-K (z-a 



(41) 



Note again that the single transverse spin asymmetry is due to the oddcron- mediated interactions with the target. 



D. STSA in Gluon and Photon Production 



Having laid out the methodology in Sections III Bl and III C[ we can now perform similar calculations of STSA for 
the cases of gluon and photon production. 

Wc begin with the gluon production. The gluon production diagrams arc shown in Fig. 1121 Since now we tag on the 
gluon, its transverse-space positions are different on both sides of the cut, now denoted z and y, while the untagged 
quark has the same transverse positions x in the amplitude and in the complex conjugate amplitude. We see that 
to obtain the gluon production cross section from the quark production expression found in the previous Section, we 
need to interchange 



z O x and y -f-» x 



(42) 



in the wave function and its complex conjugate correspondingly. In addition, since we are interested in the differential 
cross section per unit gluon rapidity yo, we use 



dye 



a 



1 



dy q 



(43) 



(with a still the fraction of the incoming quark's longitudinal momentum carried by the final state quark). 

y 




+ 




+ c. c. 



FIG. 12. Diagrams contributing to the gluon / photon production cross section. 

These modifications lead to the following expressions for the polarization-dependent and unpolarized cross sections 
for gluon production 

C F 



(2^ 



J d 2 x d 2 y d 2 z e - ife -C*-«) <f> pol ( x -z,x-y,a) X ( ^ tl (x , y , z) 
dvWp = 2(2&3 / d 2 xd 2 yd 2 ze- lk(z - v) <S> unp {x - z , x - y, a) 2^l m (x , y , z) 



(44a) 
(44b) 
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with $ po i and & unp still given by Eqs. (fTS"]) and (| 14|) . 

The interaction with the target for the gluon production case can be calculated along the similar lines to the above 
calculation of quark production by using Fig. \12\ yielding 

_ / ? n c ~ ~ 1 N c - - 1 Nc n f 1 f 

~ X^"" + 2Cf" ^ y ^ z ~ 2iV c C F " 2^ ^ + 2iV c C F ^ ~ 20p uy yx + 2N^ ° u 

(45) 

Note that now 

u = x + (1 — a) (z — x) (46a) 
w = x + (1 — a) (y — x) (46b) 

due to the interchanges of Eq. ([42]) carried out in Eqs. (|26|) . 

Separating the interaction into the symmetric and anti-symmetric components under the z -f> y interchange one 
obtains 

Nc/aa a A \ 1 /a a a a \ . 1 

Su y Sy x U y Oy X J ~T~ — S u x J , (47&) 



2 ^ m -(s»« + 2 J (s*y^* 6 zy d yz ) 2NC 2 "(s xz s zw 6 XZ 6 ZW ) + 2N s, 



2C F V — y y " ~-y~y-j 2N C C F 



Z C^ 1 \ / Z iV c O 

k Su y Oy x ~\~ U y Sy X I ~^ ~c^~\J X / ■) (47b) 



2C F V — y y - ~"y~»"j 2N C C F 

where we have used the fact that O yz = —O zy which follows from the definition in Eq. (|29c|) 
Finally, in the large- N c limit Eqs. (|4"T)) simplify to 



synim 



j(G) 
anti 



Bit W ~l~ (S Z y ) S X Z S Z W S U y Sy X [O Z y ) ~|~ X Z O Z W ~\~ U y Oy X (48a) 

large-A^ 

i [Ou w S X Z O Z w X 2 S Z 11} Su y Oy X O u y Sy a;] • (48b) 

large- -ZV C 

Eqs. (|4"T)) and (|44p . when used in Eq. ([!]), give an expression for the gluon STSA in the CGC formalism. This is 
another main result of this work. 

Constructing the cross sections for prompt photon production out of the gluon production cross sections we have just 
derived is straightforward. One has to drop all color factors in the light-cone wave functions, replace a s — > a. em Z 2 
with Z f the electric charge of a quark with flavor / in units of the electron charge, and recalculate the interaction 
with the target remembering that the photon, in this lowest order in olem approximation does not interact. One 
obtains the polarization-dependent and unpolarized cross sections for photon production 

d(A<r«) = JL J d 2 xd 2 yd 2 ze- ik < z -^% ol {x-z, X -y,a)X^ u (x,y,z) (49a) 

d °ui = ^3 / d 2 xd 2 yd 2 ze- ik <*-^ <S> unp ( x -z,x-y,a) l^ mm (x , y , z) , (49b) 

where & P oi and $ UTtp are given by Eqs. (|15j) and (|14j) with the a s — > olem Z 2 replacement. 
The interaction with the target is calculated to be 

lW=l+fl„-5„-fl„ (50) 

with the symmetric and anti-symmetric under z <-> y parts 

-^syrnm- ^ ^uw S xw S ux (51a) 
-^anti = * [Ouw ~ O xw — O ux ] . (5Vo) 

Eqs. ([5"Tj) and (|4^|) along with Eq. (TT]) give us the prompt photon STSA. This is the third and final main formal result 
of this work. Note that below we will show that Eq. (|49ap leads to d(Aa^') = for any target, which implies zero 
STSA for photons in our mechanism. 
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We have constructed general expressions for STSA generated by quark, gluon, and photon production in + A 
collisions. Knowing the light-cone wave functions squared (O, (|T5|) and the interactions for the 3 channels ([50)1 
(along with (j4lj)). (|47|) . (f5~i"T) one can make explicit predictions for the corresponding asymmetries. In general terms, 
we have shown that in this formalism the asymmetry is generated by coupling of the spin-dependent part of the wave 
function to the Odderon interaction with the target. 



III. EVALUATIONS AND ESTIMATES OF THE ASYMMETRY 



Unfortunately, Eqs. (|32|) . (|44|) . and (|49| are too complicated to be integrated out analytically in the general case. 
In this Section, in order to understand the qualitative behavior of our results, we evaluate the integrals analytically, 
taking the interaction with the target in the quasi-classical Glauber-Mueller approximation. In such a quasi-classical 
limit, the real part of the S matrix (|29b[) is [75| 



Jx y 



exp 



v\ 2 Ql 



In- 



1/1 A 



(52) 



where the quark saturation scale scale Q 2 (b) is defined in terms of the nuclear profile function (transverse nuclear 
density) T(b) as 



1(b) 



Ac 



T(b) 



(53) 



and A is a non-perturbative IR cutoff. 

In the same quasi-classical approximation the odderon amplitude is [7l| 



Ox y — 



with the constant (721473 ^ 



c a 3 s In 3 



\x — r 



exp 



CO 



(A c 2 -4) (A c 2 -1) 
12 A, 3 



In. 



y A 



(54) 



(55) 



The logarithm cubed in Eq. (J54J) arises due to the triple gluon exchange between the dipole and some quark in the 
target nucleus located at transverse position r. Angle brackets in Eq. (j5~4")l denote the averaging over positions of the 
quark in the nuclear wave function, along with the summation over all the nuclcons in the nucleus that may contain 
this quark. This averaging is carried out in Appendix |Al yielding 



O* 



3 3 7T 2 

-coa s — \x-y\ exp 



I a - y\ 



in 



i 



x-y\A 



{x-y)-VT[- 



(56) 



For simplicity we will also work in the large- A c limit for the light-cone wave function. Just like before, we mainly 
concentrate on the quark production case in Eq. (|32p : STSA in the gluon production channel can be evaluated along 
similar lines. We will also consider STSA for the prompt photon production. 



A. Single Transverse Spin Asymmetry in Quark Production 



1. Spin- Dependent Quark Production Cross Section 



First let us evaluate the numerator of the STSA in Eq. ([I), which, in the quark production case, is given by 
Eq. (|32a[) . (For simplicity we assume that a < 1 which allows us to drop the contribution from Eq. (|4T)|) .) Working 



3 Note that the sign is different from that in WX \fi : the sign in Eq. (1 55 I t arises when using a consistent convention for the sign of the 
coupling g both in the Wilson lines and in the classical gluon field of the target. (Our sign convention is to have +i g for the quark-gluon 
vertex, resulting in +i g in the Wilson lines 1161 .^ While the physical conclusions reached in [7ll474H are independent of the sign of the 
odderon amplitude, the direction of the asymmetry in question explicitly depends on the sign of O xy . 
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in the large- N c limit for the light-cone wave function we substitute the interaction from Eq. (|31b|) into Eq. (|32a[) to 
obtain 

d(Aa (9) ) = i [ d 2 x d 2 y d 2 z e -^ z ^ $ , (z - x , y - x,a) [O zy + O uw - O zx S xw - O ux S xy 

z [/tt )° i — a J 

Sz X Ox W S u x X y\ . (^^) 

Our goal now is to evaluate this expression using the S'-matrix from Eq. (|52[) and the odderon amplitude (|5GI) . 

The interaction with the target in Eq. (|5T[) is non-linear. It is tempting to try to simplify the problem by neglecting 
all the multiple rescattering saturation effects. In such a linearized approximation Eq. (|57|) reduces to 

d(A<r(«)) Hn = *- ? ^3-^- / d 2 xd 2 yd 2 ze- ih < z -y^ pol (z-x,y- X ,a) 
I l — a J 

x [o z y + o u w — o zx — o u x — o x w — o x y] (58) 

where 



C \x - y\ 2 (x-y)-VT[ ) (59) 



is the linear part of the averaged odderon amplitude (|5"f)|) . However, one can easily show that the cross section in 
Eq. (|58|) is in fact zero, i.e., that 



d(Aa^) Hn = 0. (60) 
We illustrate this by considering the o zy term in Eq. (|58[) . Defining new transverse vectors 

z = z-x, y = y-x, (61) 
we rewrite the o z y contribution to the cross section in Eq. (|58| as 

l 2^T^ J d 2 yd 2 2e- ik <*-^<f> pol (z,y,a) J d 2 x o z+x , v+x . (62) 
This expression is zero since 



J d 2 x o z+XtV+x = (63) 



due to the fact that the odderon amplitude (|29c[) (and, therefore, the linearized odderon amplitude ([59)) ) is an anti- 
symmetric function of its transverse coordinate arguments, 

O xy = -O yx . (64) 

The argument goes as follows. Employing Eq. (|64j) and shifting the integration variables we write 

f(v) = J d2x Ox. x+y = - J d 2 x O x x - y = -f(-y). (65) 

Since f(y) depends only on one vector y and is a scalar under the rotations in the transverse plane, it is a function 
of y 2 only, and can satisfy ([65]) (i.e., can be an odd function of y) only if f(y) = 0. This demonstrates that 

d 2 xO XtX+y = 0. (66) 

Similar arguments can be carried out for other terms in Eq. (|58j) . leading in the end to Eq. (|60|) . We arrive at an 
important conclusion: STSA cannot result from the interaction with the target mediated by the odderon exchange 
alone. Neglecting the interactions contained in the dipole ^-matrices in Eq. ([57]) would lead to zero transverse spin 
asymmetry. This is an important observation elucidating the nature of our result (|32j) and the corresponding STSA: 
in order to generate a non-zero STSA the interaction with the target has to contain both the C-odd and C-even 
contributions! 



18 



Returning to the general case of Eq. ((57)) we see that the argument we have just presented demonstrates that the 
O zy and O uw terms are zero in the general case as well, since they are not multiplied by the S'-matrices. Dropping 
these terms yields 

d(Aa^) = -i _2L- f d 2 xd 2 yd 2 ze- lk{z -y^ <P pol (z -x,y-x,a) 

2 (2Tr) fi 1 — a J 

X [O z x S x w ~|- O u x S x y ~t~ O x w S z x ~t~ O x y S u x ] ■ (^^) 

To evaluate Eq. (|67|) let us first study its large-fer asymptotics. Since m < m and the quark mass m is at most the 
constituent quark mass of about 300 MeV (we assume light quark flavors), we have kx 3> Q s 3> rh- Changing the 
coordinates using Eq. (|61[) reduces it to 
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(68) 



For each term in the square brackets of Eq. (|68|) the integrals over z and y factorize: taking the large-fey limit in each 
of them separately, we see that the large-fey asymptotics corresponds to small Zt and yr- We thus need to expand 
the interaction with the target term in the square brackets of Eq. (|68p to the lowest non-trivial order in zr and yx- 
Note that above we have seen that if we keep the dipole iS-matrices at the lowest order in the dipole size, 5=1, 
then the spin-dependent cross section would be zero. We thus use Eqs. (|52|) and (|56|) to expand the 5-matrices to 
the next-to-lowest order, while keeping the oddcron amplitudes at the lowest order given by Eq. (|59[) . Performing the 
expansion, substituting the wave function squared from Eq. (|15[) (also expanded to the lowest non-trivial order in 2t 
and into Eq. (|fJ5)l . and employing Eq. (|53p we obtain 
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(69) 



Since zt and j/t are small, one may think of neglecting them compared to x in the arguments of T's in Eq. (|69[) . 
However, this would lead to a zero answer after integration over x. The reason for this conclusion is that any 
unpolarizcd target, after averaging over many events, is rotationally symmetric in the transverse plane. This implies 
that VT(s) — VT(ir) = xT'(xt) where x is a unit vector in the direction of x and T'(xt) = dT(xr)/dxT- 
Integrating x over the angles of x would give zero. 

Instead of neglecting zt and yr, we shift x x — z/2 in the first term in the square brackets of Eq. (|69|) and 
expand T along the lines of Eq. (|A7|) . and perform similar operations to the other terms in the brackets obtaining 
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where we have dropped the tildes over y and 2, since now it would not cause confusion. 

Using VT(x) = VT(it) = xT'(xt) and integrating over the angles of x reduces Eq. ([70)) to 



(70) 
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Integrating over y and 2 in Eq. ((71)) and discarding delta-functions of k (since fey 7^ 0) yields 



a zt A 



(71) 
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(72) 



We see that the polarized spectrum falls off rather steeply with fcr, scaling as l/fc T . This indicates that in the standard 
collinear factorization framework our STSA generating mechanism originates in some higher-twist operator. 
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Another important qualitative feature one can see in Eq. ([72]) is that the spin-dependent cross section falls off with 
decreasing longitudinal momentum fraction a, which implies that the corresponding STSA decreases with decreasing 
Fcynman-a; of the projectile, in qualitative agreement with the experimental data. 

To improve on Eq. (f72"]) let us find the spin-dependent differential cross section d(Aa^) for lower kr, closer to the 
saturation scale. To be more specific let us relax the kx 2> Q s restriction and consider a broader region of kx < Qs 
and kx ^ Q s , but still with kr S> fh. For such not very large &t we can neglect the logarithms in the exponents of 
Eqs. (|52[) and (|56p as slowly varying functions compared to powers they multiply [56l. l76l l87j. writing 
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(74) 



Substituting Eqs. ([73]) and {72]) into Eq. 
yields 



, expanding the polarized wave function squared, and dropping the tildes 
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(75) 



Similar to the large-fcr asymptotics, we shift x ^ x — z/2 in the first term in the square brackets of Eq. (|T5[) and 
expand the resulting exponential with the help of Eq. (1531) as 



4 4 Q 2 s ( X )-\ a? 4 Q^+SLVf^) 
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(76) 



The 1 in the square brackets of Eq. (|76l) does not contribute as its contribution vanishes after integration over the 
angles of x in Eq. (|75[) . leaving only the second term to contribute. Performing similar expansions in the other terms 
in the square brackets of Eq. (|75|) we obtain 
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(77) 



where we have also dropped lnl/mi/T and lnl/m^T, since, with our precision, similar logarithms were neglected in 
Eqs. (fT3"|) and ([72)) above as slowly varying functions of their arguments 0- Again, integrating over the angles of x 
yields 
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Integrating over y and z we get 
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(79) 



This is the final expression for the STSA-generating cross section for quark production. Note again that k 2 is the 
y-component of the quark's transverse momentum k = (k 1 , k 2 ). Let us point out a few of the important features of 
Eq. ([79]). First of all, we see that similar to Eq. ([72]) decreases with decreasing a for small a, now due to the factor 
of 1/a 2 in the exponent. We also see that for kr — > the cross section d(Aa^) also goes to zero. We also note that 
the cross section ([79]) is not a monotonic function of fcy. In particular, for positive k 2 it starts out negative at small 
kr, becoming positive for kr > Q s a y/l + a 2 / (1 — a 2 ), in agreement with the large-fcr asymptotics of Eq. (|72[) . 



4 We have done the calculation without neglecting those logarithms: the resulting changes were mainly of quantitative nature, while the 
obtained expression was significantly more complicated than Eq. 1791 . Since both the expressions with and without the logarithms are 
approximate, we decided to only show the latter in this work due to its relative compactness. 
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2. The Unpolarized Cross-Section 

The real hadronic STSA in Eq. ([1]) contains contributions from both quark and gluon production cross sections 
(|32p and (|44[) in the numerator and in the denominator, convoluted with the fragmentation functions for the quarks 
and gluons decaying into a particular hadron species as well as the transversity distribution of polarized quarks. This 
is what needs to be done to have a real comparison of the data with our theoretical results. While such comparison is 
beyond the scope of this work, we would like to assess the main qualitative features of our STSA-generating mechanism 
by concentrating on quark STSA only. 

It may be tempting to consider a situation where both the numerator and the denominator of Eq. ([1} are driven by 
the quark contributions. However, the unpolarized valence quark production cross section (|32b[) is known to decrease 
with decreasing quark momentum fraction a [77l [78j , while both the unpolarized gluon and sea quark production cross 
sections grow with decreasing a in theoretical calculations [HI, |89[ . In the actual experiments the hadron multiplicity 
also increases as we move further away from the projectile in rapidity. 

Therefore, in order to get a somewhat realistic evaluation of the qualitative behavior of the obtained STSA, we 
will use the unpolarized gluon production cross section in the denominator of Eq. ([T|). While evaluation of the 
unpolarized gluon cross section (|44b|) along the same lines as were used to obtain Eq. ([75)1 is somewhat involved, we 
will approximate the result by assuming that the produced gluo n is soft (i.e., far from the projectile in rapidity), in 
which case the corresponding production cross section is [5a "l76l. \8l\ 

, lr s a s N c 7 , 9 [ 1 2 1 fi— \ ( kl \ , 4 k 2 A 2 1 1 
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3. Single Transverse Spin Asymmetry 

We now have all the essential ingredients to sketch the STSA due to quark production in the large- N c limit (for the 
wave function): we have Eqs. (|79[) and (|80p . giving the numerator and the denominator of Eq. ([1]) correspondingly. 
We thus write 
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(81) 



The xt~ and j/T-integrals in Eq. (|8ip appear to be very hard to evaluate analytically. Instead we evaluate the integrals 
numerically assuming a simple Gaussian form of the nuclear profile function, 



T{b) = \Rpe- b -l R2 



(82) 



with R the nuclear radius and p the nucleon density. Such Gaussian profiles are of course not realistic for nuclei, but 
have been successfully used to describe protons (see e.g. [90j). 

In evaluating the STSA in Eq. (|8"Tj) one has to remember that in the standard convention one has to choose k in 
the direction left of the beam, which, in our notation, means along the negative y-axis (see Fig.[T]). Hence we need to 
replace k 2 — > —kr in Eq. (I8T1) . 

To plot Eq. (|8"Tj) we will attempt to use somewhat realistic numbers, while realizing that all the theoretically- 
calculated cross sections are likely to have non-perturbativc normalization corrections, which may affect the size of 
the effect. To that end, we will use the saturation scale (cf. Eq. ([53]) ) 



Q 2 (b) = 2™ 2 K 2 T{b) 



(83) 



with the iT-factor fixed at K = 10 to make Q, rj 1 GeV, which is a realistic value for a proton. (Each T'{xt) in 
Eq. ([81) is multiplied by the same if 2 -factor, since it also arises from the saturation scale.) We put m = 300 MeV 
to mimic a constituent quark, along with p = 0.35 fm -3 for a proton of radius R = 0.878 fm, and a s — 0.3. We plot 
the resulting A$ from Eq. (f8"Tj) in Fig. [T^Jfor different values of a with the IR cutoff A = 100 MeV and cutting off 
the Xt~ and yr integrals in Eq. (fBTj) at 2.1 fm in the IR. (Note that strictly-speaking the CGC formalism employed 
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here is valid only for scattering on a nuclear target, since it resums powers of a large parameter a 2 s A 1 / 3 . However 
its applications to proton target have been successful phenomenologically in the past 
estimates here could be relevant for p' + p collisions.) 



91], giving one hope that our 




FIG. 13. Quark STSA from Eq. (|81|) for the proton target plotted as a function of kr for different values of the longitudinal 
momentum fraction a carried by the produced quark: a — 0.9 (dash-dotted curve), a = 0.7 (solid curve), a = 0.6 (dashed 
curve), and a = 0.5 (dotted curve). 



From Fig. 1131 we see that our STSA is a non-monotonic function of transverse momentum fcx, first rising and then 
falling off with kx in qualitative agreement with the data shown in the right panel of Fig. As one can clearly see 
from Eq. (|81[) the maximum of A$ at impact parameter xt in our formalism is determined (up to a constant) by 
the saturation scale, ~ Q s {xt)j such that the asymmetry integrated over all impact parameters peaks at kj- ~ Q s 
with Q s an effective averaged saturation scale. The conclusion about An peaking at kx ~ Q s was previously reached 
in [62j. Let us stress again that the STSA in our case changes sign when plotted as a function of kx or a (i.e., it has 
a "node"). 

Note that, while the magnitude of STSA plotted in Fig.[l3]can be as large as tens of percent, like the data in Fig. [21 
the momentum at which the asymmetry is non-zero appears to be much smaller in our Fig. [13] than it is in the data 
of Fig. [2] The discrepancy of the /cy-range of the data and our Fig. [T3] signals the following potential problem: the 
XT-integral in Eq. ([81]) is dominated by large xt, where Q s {xt) is small, leading to small values of kr dominating 
Aff, and potentially making the corresponding physics non-perturbative. Thus our perturbative calculation appears 
to be sensitive to the non-perturbative domain. 

To illustrate the range of spectra that can be obtained by our estimates, we rcplot An from Fig. [T5] in Fig. [TJ] 
cutting off the xt- and ?/T-mtegrals in Eq. (|81[) by 1.3 fm. In addition, we mimic the coordinate-space logarithms, like 
those that were neglected after l|75[) . by introducing a factor of hikx/m- In this plot the fc^-range of the asymmetry 
is broader than in Fig. 1131 which makes it closer to the experimental data in Fig. [5J but the height of the asymmetry 
is over an order-of-magnitude lower than the data. More work is needed to assess whether the cutoff dependence is a 
result of the approximations made, or whether it actually signals a potential breakdown of the approach indicating 
the non-perturbative nature of STSA. 

Another important observation one can make from Fig. [13] is that A$ increases with increasing a, except for very 

large values of a when it starts to decrease. The increase of A% ] with increasing a is in qualitative agreement with 
the data in the left panel of Fig. [5] where the data points increase with increasing Feynman-x. One can also see that 
while our STSA in Fig.[13]is mostly positive, the data for ir~ mesons in the left panel of Fig.[2]gives a negative STSA. 
This result can be explained in the following simple model. Imagine a constituent-quark model of a proton, with the 
spins of both up quarks aligned with the net proton spin, and the spin of the down quark pointing in the opposite 
direction. Then in our mechanism each up-quark would give a positive STSA denoted A u , while the down-quark 
would give a negative STSA of the same absolute value, denoted Ad = —A u . Imagining that after the collision the 
up-quarks fragment into 7r + 's and 7T°'s, while the down quarks fragment into 7T°'s and 7r _ 's, and taking into account 
that there are twice as many up-quarks than down-quarks in the proton, and neglecting pions resulting from the gluon 
fragmentation, we obtain A n + = A u , A n a = A u /2, and A v - — Ad/2 = —A u /2, thus obtaining a negative STSA for 



22 




FIG. 14. Same as in Fig. 1131 but with the 1.3 fm upper cutoff on the xt- and yr integrals in Eq. (|81[) and a factor of Infer/™, 
inserted. 

7r - 'slf| This naive model appears to be in a qualitative agreement (and in loose quantitative agreement) with the data 
in the left panel of Fig. [2j Inclusion of gluon fragmentation may further improve the quantitative agreement with the 
data. 

Finally, to test the dependence of our STSA in Eq. (|81"j) on the size of the target, we note that for kx ~ Q s one gets 

A%\k T ^Q s )~-^~A- 7 / 6 , (84) 

if Q 2 S ~ A 1 / 3 . This indicates a very steep falloff of STSA with the atomic number of the nuclear target. Such a 
conclusion appears to be supported by the numerical evaluation of Eq. (|81[) for several different radii of the target 

shown in Fig. [151 (Now the xt- and dt integrals are cut off at 2.4 fm.) One can see that A$ drops very rapidly 
with the size of the target. If the experimentally observed STSA in p^ + p collisions are due to our mechanism, our 
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FIG. 15. Quark STSA from Eq. (1811) plotted as a function of kr for different values of the target radius: R = 1 fm (top curve), 
R — 1.4 fm (middle curve), and R = 2 fm (bottom curve) for a = 0.7. 

prediction is then that in p^ + A collisions STSA should be much smaller than that in p^ + p. In the case of a heavy 
ion target like Au the STSA due to our mechanism is likely to be negligibly small. 

While we have demonstrated here the potential for our calculations to agree with the data, the evaluations presented 
here have to be significantly improved to reach a definitive conclusion. For instance the fc^-dependence and the overall 



5 We would like to thank Mickey Chiu for a discussion of these estimates. 
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normalization in our Eqs. (|79p and (|80[) are overly simple and is not ready to be compared to the data. The equations 
need to be corrected for the effects of DGLAP evolution, small- a; evolution and for the running of the coupling for 
a meaningful quantitative comparison with the data. Only such a phenomcnological analysis can determine whether 
our mechanism for generating STSA is dominant, or whether it is simply one of the many factors contributing to the 
asymmetry. 

B. STSA in Photon Production 

Using the methods developed in Sec. IIII Al we can now evaluate the photon STSA given by Eqs. (|49[) and (|5ip . 
Substituting Eq. (|51b|) into Eq. (|49a[) and performing the variable shift of Eq. (|6ip while keeping in mind that now it 
and w are given by Eqs. (|46|) yields 

d(Acr (7) ) = j^-^ J (Fx<Pyd?ze~ lh '( z ~ y ) $ po i(-z, -y,a) [O as+ (i-eO je , SB +(i-a i ) s - XtX+{1 _ a)y - O^i-a)*, «] 

(85) 

where we have again dropped the tildes for brevity. Using the argument of Eqs. (|65|) and (|66j) we see that each term 
in the square brackets in Eq. (|85|) is zero after the integration over x. We thus have an exact result that 

d(Aa (7) ) = (86) 

in our mechanism for generating photon STSA. Hence the photon STSA is zero, = 0, in the forward production 
region under consideration. 

IV. CONCLUSIONS 

To conclude let us summarize the main points of this work. Above we have shown how STSA can be generated 
in the CGC formalism for quark and gluon production. The results for the corresponding cross sections are given in 
Eqs. (f3"2l. (|4"U)) and (|4~4")l . The same mechanism gives zero STSA for prompt photons. 

In our case STSA is generated by both a splitting in the projectile wave function, and by the combination of the 
C-odd and C-even interactions with the target. Hence our STSA-gencrating mechanism is distinctively different from 
the Collins [13] and Sivers [H|,[l6| effects, and is more akin to (though still different from) the higher-twist mechanisms 

of 

Evaluating the quark STSA in a simplified quasi-classical model we found qualitative agreement with the data: 
quark STSA appears to be a non-monotonic function of kx, and is an increasing function of increasing xf (for most 
of the xp-range). It is perhaps encouraging that the obtained asymmetry can be of the order-of-magnitude of the 
experimental data. Further phenomcnological studies of our formulas (|32[) and (|44[) are needed to determine whether 
this qualitative agreement of our results with the data can become quantitative. 

Analyzing the general quark production formula (|27[) one can see that the contribution to STSA arises from the 
z <-> y anti-symmetric part of the integrand. In arriving at Eq. (|32a|) from Eq. (|27|) we employed the lowest-order 
(order-a s ) spin-dependent part of the light-cone wave function squared (|15p . which happens to be symmetric under 
the z -o- y interchange: hence, in our case, to obtain a contribution to the STSA the interaction with the target had 
to be z <-> y anti-symmetric. However, it is possible that higher-order corrections to the light-cone wave function 
squared would lead to a z <-> y anti-symmetric contribution. (By the 'wave function corrections' we understand all the 
initial and final state corrections with rapidities between the projectile and the particle we tag on.) In such case, the 
interaction with the target need not be z y anti-symmetric, and can be mediated by the standard C-even exchange. 
To test whether such scenario is feasible within the CGC/saturation perturbative framework one has to calculate the 
higher order corrections to the polarization-dependent light-cone wave function squared ()15|) . The corrections would 
need a generate a relative complex phase between the corrected and uncorrected wave functions, in agreement with 
the initial proposal of [2^, [H|. This can be accomplished in LCPT if the corrections lead to an intermediate state, 
for which the imaginary part of the energy denominator leads to a non- vanishing polarization-dependent contribution 
to the scattering amplitude. An example of such corrections in our case could be a modification of the amplitude 
in Fig. [5] resulting from a gluon exchange between the outgoing quark and gluon formed in the projectile splitting. 
Calculation of such diagrams appears to be rather complicated and is beyond the scope of this work. However, it 
potentially may give a contribution comparable to the STSA resulting from Eq. (|32a|) : the latter consists of the 
order-a s light-cone wave function squared, convoluted with the target interaction resumming powers of a 2 s A 1 / 3 and 
a s Y, with one extra power of a s due to the odderon exchange (|56p . Our contribution Q32a[) is, therefore, order-a^, 
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if one assumes that a 2 s A 1 / 3 ~ 1 and a s Y ~ 1, which is parametrically comparable to the C-odd order-a^ light-cone 
wave function squared, interacting with the target through a C-even order-one exchange. An explicit calculation is 
needed to explore this possibility and is left for future work. 
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Appendix A: Averaging the odderon amplitude over target wave function 



Let us construct the dipole odderon amplitude averaged over the target field. The triple gluon exchange happens 
between the dipole and a nucleon in the target, which, for simplicity we model as a valence quark in a bag. (The overall 
factor in front of the averaged odderon amplitude should indeed depend on the details of the averaging; however, we 
believe the coordinate-space dependence would remain the same for other models of the nuclear wave function.) The 
target averaging then consists of averaging over the positions of the quark in the nucleon and over the positions of 
nucleons in the nucleus, along with summation over all nucleons. Assuming, again for simplicity, that the quark has 
equal probability to be anywhere inside the nucleon in the transverse plane (a cylindrical "nucleon" approximation), 
we write for the averaged odderon amplitude 
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2 J ~ |as-y|Aj 

Here b is the position of the center of a nucleon in the transverse plane with respect to the center of the nucleus, r is 
the position of the valence quark in the nucleon, and a is the radius of the nucleon, as illustrated in Fig. 1161 The two 
thcta- functions in Eq. (|A1[) insure that the valence quark and the center of the a:,y-dipole are both located inside 
the nucleon in the transverse plane. In our simple model of the collision the dipole has to hit the nucleon directly in 




nucleus 

FIG. 16. The geometry of the dipole-nucleus scattering as employed in Eq. (jAljl . 



order to be able to interact with the quarks inside of it. Since the dipole x, y is perturbativcly small, we enforce this 
condition by demanding that only the center of the dipole is inside the nucleon's transverse extent. 
To integrate over r in Eq. (|Aip we notice that 

(A2) 
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if the integration carries over the whole transverse plane. Using this result wc write 



d 2 r ln d 



x 



\y-b 



(a — r) = — / d r In' 



\x — b — r | 
\y-b-r\ 



(r-a). 



(A3) 



25 



To approximate the integral on the right- hand- side of Eq. (|A3[) we expand its integrand in powers of \x — b\/r and 
\y — b\/r to the first non-trivial (after integration) order, thus obtaining 



d 2 r In 3 



\x — b — r\\ „. 37r . l9 , , . 

y -6-r j e{a-r)^— 2 \ X -yf (x y) ■ (x + y 2 b). 



Substituting Eq. (|A4[) back into Eq. (|A1|) and defining a new integration variable 

a; + y 



yields 
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(A6) 

In principle this result is as far as one can simplify O xy without the explicit knowledge of the nuclear profile 
function T(6). To obtain a closed expression for the STSA in the text we expand 



T 



y 



T 



y 



(A7) 



with V the transverse gradient operator. Such expansion is potentially dangerous near the edge of the nucleus profile, 
where the derivatives may get large. For instance, for a solid-sphere model of the nucleus the nuclear profile function 
is T(b) = p 2 V-R 2 — b 2 with p the nucleon density and R the nuclear radius; the derivatives of such T(6) near b = R 
are divergent. Using the realistic Woods-Saxon profile would make the derivatives finite, but they would still be large. 
Thus we will proceed by using the expansion (|A7[) as a way to simplify the expression, keeping in mind that in the 
cases where this expansion breaks down one has to return back to Eq. (|A6|) . 

Substituting Eq. (|A7[) into Eq. (| A6|) and integrating over 6 yields (for the first non-trivial term after integration) 
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(A8) 



This is our final expression for the target-averaged odderon amplitude. Note an interesting feature of Eq. (|A8|I : the 
non-zero contribution to the odderon amplitude in the transverse coordinate space arises from the gradient of the 
nuclear profile function. The odderon interaction with the target is thus only possible if the target has a non-uniform 
profile in the transverse space. This is in stark contrast to the C-even exchanges, which are non-zero even for the 
Bjorkcn model of a nucleus of infinite transverse extent with constant density in the transverse plane. 
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